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Quantum corrections coming from massless fields conformally coupled with gravity are studied, in order 
to see if they can lead to avoidance of the annoying Big Rip singularity which shows up in a flat Friedmann- 
Robertson- Walker universe filled with dark energy and modeled by a scalar phantom field. The dynamics of 
the model are discussed for all values of the two parameters, named a > and /3 < 0, corresponding to 
the regularization process. The new results are compared with the ones obtained in 1 1 1 previously, where dark 
energy was modeled by means of a phantom fluid with equation of state P = top, with uj < — 1. 
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I. INTRODUCTION 

Recent observations of distant type la supernovae, baryonic acoustic oscillations (BAO), anisotropies of the cosmic 
microwave background radiation (CMB), and some other, confirm that our universe expands in an accelerated way J2JE1- 
In fact, it seems that it is at present in a dark energy phase 0J. A proposal to explain this situation is to assume that 
the energy density of our universe is dominated by a phantom scalar field, that is, a model where the energy density 
and pressure are p — —\(j) 2 + V(<fi) and P — —\<\> 2 — V((f), being <f> the scalar phantom field. In this case, future 
singularities are bound to appear in a finite time [5|. These singularities are undoubtedly there in the classical situation 
when no quantum effects are taken into account, but it seems feasible that, near the singularities, where the curvature 
has very high values, quantum effects could have the power to drastically modify the behavior of the universe, yielding a 
milder singularity or maybe even a non-singular model. 

In this paper we extend the study carried out in [1 J to a dark energy universe modeled by a scalar phantom field. In fact 
we will consider exponential potentials which give rise to a Big Rip singulary, and introduce quantum corrections in order 
to avoid these late time singularities. Specifically, we shall consider the quantum effects due to massless, conformally 
coupled fields. This is a special, workable case where the quantum vacuum stress tensor — which depends on two regu- 
larization parameters, here called a > and /3 < — and the semiclassical Friedmann equation, can be both calculated 
explicitly. 

We will show, analytically and numerically, that quantum effects drastically modify the Big Rip singularity, rendering 
it of type III or turning it into a singularity in the contracting phase. In the first case (a type III singularity) the Hubble 
parameter does not diverge, but the energy density does tend towards infinity. In the other case (a singularity in the 
contracting phase) the Hubble parameter becomes finite and negative, and the energy density diverges towards minus 
infinity. What is important to note is that, in both cases, the Hubble parameter remains finite. 

The paper is organized as follows. In the next Section, using the mathematical theory of dynamical systems, we study 
some phantom fields driving the universe to a Big Rip singularity. In Sect. Ill we introduce the quantum corrections due 
to a massless conformally coupled field, and we perform an analytic study of the semiclassical Friedmann equations. In 
Sect. IV a numerical analysis is carried out to check to good approximation the analytic results obtained in the previous 
Section. In Sect. V we analyze the problem in the context of loop quantum cosmology, where it has been stated that 
quantum corrections do completely avoid the Big Rip singularity. We will see that the way to obtain these conclusion is 
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in doubt, because they have been got in some places from an incorrectly modified Friedmann equation. In last Section we 
compare the results obtained for a phantom field with those that were derived for a phantom fluid model. The units to be 
used in the paper are: c = h = M p = 1, where M p is the reduced Planck mass. 



II. DARK ENERGY MODELED BY A PHANTOM FIELD 



A phantom field (ft is modeled by an energy density of the kind p — — ^(ft 2 + V((ft) and a pressure P 
For this field the Friedmann and conservation equations are: 

' H 2 = \(-\j> 2 + V(cft) 



(1) 







-<f>-3H<j> 



dV 

d<t> 



From this system one deduces that p = 3H<ft, and thus H = ^(ft 2 > 0, which means that ~ = H + H 2 > 0, that is, in 
this model the universe is expanding in an accelerating way. 

To analyze the dynamics of the system we start considering a power-law potential, i.e., V(<f>) = X(f> 2n with A > 0. The 
field equation can be written as follows 



d_ 



1 



+ V(4>) = -3H<p, 



(2) 



where V((ft) — —V{<j>). This is a dissipative system, and the slow-roll conditions (V'/V) <C 1 and |V"/V| <C 1 are 
satisfied when \(ft\ ^> n. Then, due to the attractor nature of the slow-roll regime, at late time, the solutions have the same 
behavior as the slow-roll solution, which satisfies the system 



H 2 



\v(<t>) 



= -3H<j) 



dV 



(3) 



Since the dynamics of the system decouples for <j) > and (ft < 0, we only consider the domain (ft > 0, where the field 
obeys the equation (ft x ' 



2n\l ^<fi™ r , which solution is 



2n(n-2)J§(t s -t) 



<ftsr{t) = (ft sr (t )e 4 



for n > 2 
for n = 2 

(ft S r(t) = <ft S r(to)+2JUt-t ) for n = l. 



f(*-*o) 



(4) 



Evaluating H sr (t) = y ^(ft™ r (t), one can see that the Big Rip singularity appears for n > 2. 

The following remark is in order. To prove, in a rigorous way, the attractor nature of the slow roll solution, we may use 
the variables x = 4= and y = ^jt^ (Q)- Then, the dynamical equation is 



dy 
dx 



= -3(1 + y 2 



1 - 



V 



V&yV 



= -3(l + y 2 ) 



2n 
6yx 



(5) 



and it is easy to verify that, for large values of x, this is the 



The slow roll solution is the curve -4*- = 0, i.e., y 
leading term of the solution. Then, since g < above this curve, and t| > below it, this definitely proves that the 
slow roll solution is an attractor at late times. 

As a different specific example, we choose V((f>) — Vbe~ 2< ^° (being Vq and <ft Q two constant parameters). Then, with 
the change of function (ft = (fto ln(ift) (now ift belongs in the domain (0, oo)), the system becomes 



H 2 = 



si 2 (" 



•3. ,/,2 



(6) 



= -M^-*P 2 )-3H(ft W 



, Vn 
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In the expanding phase (H > 0), this system can be written as follows: 



and dividing this equation by <f> — i.e., using the variable <f> as a time — one obtains the equation 



This is an autonomous first order differential equation, therefore, it can be completely studied just through the sign of the 

/2Vp V2V ( 

00 ' 00 



function F. From the Friedmann equation, one can see that the domain of F is the interval [— v T i7 ° , v ?^° ]. The zeros of 



F are the points ± and —^^--j===. F has a vertical asymptotic at zero. Finally, F is positive in the interval 
( - 7 \ , . 0) and negative in ( - ^ , - 4^ . \ , ) U (0, ^ ) . 

This all means that, using <p as time, the critical points ^^J" and — ^^ 3 ^ 2 are repellers and the critical points 
and — are attractors (see Pict. 1). However, in the domain ip < 0, when the time <j> increases, the cosmic time 



t decreases, and vice versa, what means that, in terms of the cosmic time, the critical point — v , , = is a global 

attractor while the other two critical points ± are repellers (see Pict. 2). 
In terms of the field 0, the critical points obtained above are: 




t.-r 



0(t) = o lnf^^-); ff(t)=0, (9) 



2V 



where i s is an arbitrary constant. Then, since — , 1 3 is a global attractor, it follows that, except for the solutions 

ip = ± V ^J" , all the other have a Big Rip singularity. 



time 4> 



00 



00 



Picture 1 



time t 



Picture 2 
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Figure 1: Phase space portrait for cj>o = Vq = 1. Blue curves are the orbits of the system and the red ones are the level curves of H. 
The line ip = — % is the global attractor orbit. 



III. QUANTUM CORRECTIONS 



In this section we will study in detail what is the change in the dynamics in the model V(tf>) — Voe™ 2< ^° when one 
takes into account quantum effects. It is well-known that for a massless, conformally coupled field, the anomalous trace 
is given by 0][6) 



T vac = aUR- P -G, 

with R = 6 (^H + 2H 2 ^j the scalar curvature and G = 24H 2 (^H + H 2 ^j the Gauss-Bonnet curvature invariant. 
The coefficients, a and /?, coming from dimensional regularization are |£S) 

1 



(10) 



a 



/3 



-^ 0l ^( N o + ^ 1/2 + 12N 1 )>0, 
— ^(N a + UN 1/2 + 62^) < 0, 



2880tt 2 



(11) 



being No the number of scalar fields, N1/2 the number of four-component neutrinos, and N\ the number of electromag- 
netic fields, respectively. 

In terms of the Hubble parameter, Eq. ( fT0] > is [9 | 



T vac = 6a(H + 12H Z H + 7HH + AH 2 ) - 12/3(H 4 + H H). 

With the trace anomaly being T vac — p vac — 2>P vac and, inserting ( p"2] ) into the conservation equation, p v 
Pvac) = 0, the modified energy density reads 



(12) 



3H(p vac 



p vac = 6a(3H 2 H + HH - -H 2 ) - 3/3H 4 



and the semiclassical Friedmann equation becomes 



TT 2 P ' Pvac 
ti = , 



(13) 



(14) 



with p 



^ 2 + V e- 
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Using the dimensionless variables t — H + t, H = H/H + , Y = H jH\, t/> = ip and p — ^jpr, with H + — yf—l/P, 
the semiclassical Friedmann equation and the conservation equation can be written as an autonomous system: 

6aH 2 Y + aY 2 + BH 4 ) , 

(15) 



( H' 


= Y, 




Y' 


1 

2aH 


(-(3H 2 + pp- 




= f, 






v- 


-3Jf£-2-*%, 



where ' denotes derivative with respect to the time t, and we have defined the new parameters Vq = Vq/H? and <fio — <j>Q. 

What we see at first sight from this system is that it does not have any critical point. It is also easy to show that the 
energy density p evolves in accordance with the equation p' — Hcf) 2 ^, which means that the energy density increases in 
the expanding phase, and decreases in the contracting one. 

Now, we look for singular solutions of the system with the following behavior near to the singularity ifTOlfTTIl 

$(t)=A(t s -t) + B(t s -t) 2 + 0((t s -i) 3 ); H(i) = H + 5H{i), (16) 
where A, B and Hq are some constants. Inserting these functions in the conservation equation, one obtains 

A=^P>; B = -%H Q A, (17) 
and inserting them in the semiclassical Friedmann equation and retaining the leading terms, one gets 

5H"{t) = J^JL => SH(t) = - t) In ((i s - t)/T) , (18) 

la t — t s Za 

where T is an integration constant. 

What we observe here is that, when we introduce quantum corrections, the Big Rip singularity, for Ho > 0, is trans- 
formed into a type III singularity, because as i — > i s one has H — > Hq, p^co and \p\ —> oo. And, when Hq < 0, one 
gets H — > H (contracting phase), p — > — oo and \P\ — > oo. 

In order to qualitatively study the system it is quite convenient, as in HUE), to perform the variable change p = y/\H\. 
After what, the semiclassical Friedmann and conservation equations become 

-2 

p" = -d p W(p,p)-3ep 2 (p') 2 , jj = H$L, (19) 
where W[p, p) = £ (f{l - \p 4 ) + £), and e = sign (H). 



For positive values of p, the potential W (Fig. 3 of Ref . 11131 ). has a unique zero, at po = (3/2) ' ( 1 + y 1 + |p 

and two critical points, at p± = ( 1± v / i-4 g^ ^ ^ xhus, for p > 1/4 there are no critical points, being the 

potential strictly increasing, from — oo to oo. For p < 1/4, the potential satisfies W(0) = — oo and W(oo) = oo, and 
exhibits a relative maximum, at p_, and a relative minimum, at p + (a hollow one). For very small values of p, at p- one 
has H 2 = p, that is, the system is close to the Friedmann phase and, at p+, one has \H\ = 1, that is, the system is close to 
the de Sitter phase. On the other hand, for negative values of p, the potential only has a critical point at p + , and satisfies 

W{0) = W{oo) = oo. 

Now, assume that, initially, the system has an energy density which is positive, and that it is in the expanding phase 
(what does happen nowadays). Then, since in the expanding phase the energy density increases, this means that the slope 
of the potential is more steep and thus the system can evolve to the contracting phase. When it enters that phase, the 
energy density decreases and even it could be negative; if so, the system is confined in the decreasing phase because the 
potential satisfies W(Q) = W(oo) = oo. 

What is important to stress here is that the system cannot remain all the time in the expanding phase due to the form of 
the potential, and also the fact that the energy density is increasing, in this phase. Three different situations may occur: 



1. The system may develop a singularity in a finite time (type III singularity). This comes from Eq. (18). 
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2. The system may enter in the decreasing phase and the energy density becomes negative, and then the system cannot 
abandon this phase. In this situation the energy density could by — oo in a finite time. 

3. The system may bounce infinitely many times (an oscillating universe). 

This is what one can say by analytically studying the system. What we will do in next section is to perform a corresponding 
numerical study, which will to show that only the first two situations are actually possible. 



IV. NUMERICAL ANALYSIS 



In this section we numerically integrate the system ( |15) , assuming that initially the system is in the Friedmann phase, 
that is, at time t = the variables (H (0) ,Y (0) , ^(O) , <p(0)) satisfy the constrains 



H 2 (0) = l( VUh - 



Vh 



V> 2 (o) 



F(0) 



1^(0) 
2 2 (O)' 



(20) 



This means that the initial conditions depend on two variables. Next, to perform our calculations we choose as variables 
(-0(0), <p(0)), and also the following values for the parameters: Vq = 24 and (j>a — 4/V3. Note that, from Eq. (201, with 
our choose of parameters the variable (p(0) belongs to the interval [—3, 3] and V>(0) belongs to (0, oo). 

In Fig. 2 we plot three simulations for different values of (3/a (the system ( 15 i depends on this quotient), the first being 
for fj/ot = —0.5, the second for (3/a = —1, and the last one for /3/a = —10. The blue color means initial conditions 
which drive to a singularity of the form given by Eqs. ( fTT) and ( |T8] > in the decreasing phase, that is, the Hubble parameter 
is negative and the energy density diverges to minus infinity. On the other hand, the red color means initial conditions 
which drive to a type III singularity. 



In Fig. 3 we have integrated the system ( 15 1 for (3/a = —10, and we show the evolution of the Hubble parameter and 
of the energy density. In the first two plots the initial conditions are taken in the blue region of Fig. 1, given a universe 
evolving, at late times, in the decreasing phase with an energy density which diverges at late times. The last two plots 
correspond to initial conditions taken in the red region of Fig. 1, and they show a type III singularity. 






Figure 2: Three different simulations, for the values /3/a = —0.5, /3/a — —1, and /3/a — —10, respectively. Red points mean initial 
conditions driving to a type III singularity. Blue points, initial conditions driving to a singularity in the decreasing phase. 
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t t 

Figure 3: Plots of H and p for /3/a — — 1. The first two have the same initial condition, corresponding to a blue point in Fig. 1. The 
last two plots have the same initial condition, that corresponds to a red point in Fig. 1. 



V. PHANTOM FIELDS IN LOOP QUANTUM COSMOLOGY 

For the flat FRW spacetime, Einstein's theory is obtained from the Lagrangian C = ^Ra 3 + C ma tter, where a 
denotes the scalar factor and L maU er — a 3 P = fl3 ( — \4> 2 — V{4>))- This Lagrangian can be written as follows 



C = 3 ( - 1 — a 2 a\ + a 3 P, what means that the same theory is obtained avoiding the total derivative, which gives the 



Lagrangian Ce = — 3d 2 a + a 3 P. The conjugate momentum of the scale factor is then given by p = ^fi- = —6aa, and 
thus the Hamiltonian is 

dP ■ v 2 

H E = dp + a 3 ^(f>- C E = —~ + a 3 p = -3H 2 a 3 +a 3 p. (21) 
dip 12a 

On the other hand, in loop cosmology the following effective Hamiltonian, which captures the underlying loop quantum 
dynamics, is considered Il4lfl6l 



sin 2 (A/3) 

iLlqc = -3V ^ 2 \-Vp, (22) 



where 7 is the Barbero-Immirzi parameter and A is a parameter with dimensions of length, which is determined by 
invoking the quantum nature of the geometry, that is, through identification of its square with the minimum eigenvalue of 

the area operator in LQG, which gives as a result A = y ^7 (see lfl6l ). Here V is the physical volume V = a 3 and /3 is 
canonically conjugated to V, and satisfies {/3, V} = \, where {, } is the Poisson bracket. 

The Hamiltonian constraint is then given by s '" 2 ^f - = §, and the Hamiltonian equation yields the identity: 

V) m> u 1 ldU L QC M „ sin(2A/3) 1 . . 

V = {V,H L qc} = d p ^H= ^ A <^0 = — arcsm(2A7 J ff). (23) 
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Writing this last equation as H 2 = hl " 2 ^f (1 — sin 2 (A/3)), and using the Hamiltonian constraint Hlqc = 
s ™2^f' ) = 3 , one gets the following modified Friedmann equation in loop quantum cosmology 



Pc\2 
2 / 

3\' Pc) ' ' Pc/12 ' pl/A 



^4fl-^W4 + ^ = l, (24) 



being p c = -^r^z- This equation, together with the conservation equation —(f) — 3H(f) + ^ = 0, determine the dynamics 
of the universe in loop cosmology. 

From the equation of the ellipse, one can easily check that the Hubble parameter belongs to the interval [— p c /12, p c /12], 
and the energy density, p, to [0, p c ], what means that there is not Big Rip. In fact, an exhaustive study of the potential 
V = y o e" 20/0 ° was performed in ETJIISl. 

But here a problem appears. It is well-known that the current cosmological theories are built from two invariants, the 

scalar curvature R = 6 (h + 2H 2 ^j and the Gauss-Bonnet curvature invariant G = 2AH 2 + H 2 ^j. For example, in 

Gauss-Bonnet gravity lfl9l the Lagrangian Cmg — a 3 f(R, G) + a 3 P is used, and semiclassical gravity, when one takes 
into account the quantum effects due a massless conformally coupled field (see for instance HI), is based in the trace 
anomaly T vac = a\3R — % G (being a > and < two renormalization coefficients). However, from the Legendre 
transformation 

U LQC = --VP + V^j>- C lqc (25) 
7 dcj) 

one gets, in terms of the standard variables, the following Lagrangian 

3a 3 H . „ s 3a 3 



-LQC 



, arcsin(2A 7 i?) + — %— ( 1 - y/l - Aj 2 X 2 H 2 ) + a 3 P, (26) 
7A 27^A^ V / 



which is not invariant. This is in disagreement with one of the main principles of general relativity. 

From these observations, one can conclude that the modified Friedmann equation (f24]i does not stand in this form, 



because it has been obtained assuming that (22i is the Hamiltonian of the system which is in contradiction with the 
invariance of General Relativity. Finally, we also conclude that the results obtained from this modified Friedmann equation 
need a deep revision (for more details, see [20 1). 



VI. DISCUSSION AND COMPARISON WITH THE PHANTOM FLUID MODEL 

In |Q~) we have studied in detail the case of a phantom fluid modeled by the EoS P = cup with w < — 1. There, we have 
shown that, in the case — 1 < ^ < 0, there exists a one parameter family of solutions which evolves into the contracting 
Friedmann phase at late times and only a particular solution asymptotically converging towards the de contracting de 
Sitter universe. All the other solutions enter into the contracting phase and become singular at finite time, satisfying 
lim t _,. t3 H(t) = —00 and lim t „> ts p(t) — 0. On the other hand, we have shown in [TJ, that for — 1 > ^ almost all 
solutions describe a universe bouncing infinitely many times (an oscillating universe). 

In the present paper, by studying a phantom field we have shown, both analytically and numerically, that all solutions 
are singular. Some of them display Type III singularities and the other ones are singular in the contracting phase, satisfying 
H(t) -> H(t s ) < 0, p(t) -> -00 and \P(t)\ -> 00, when t -> t s . 

The difference comes from the fact that, for a phantom fluid, when one considers the dynamics in M 3 using the coordin- 
ates (H, H, p), the manifold p — is invariant. More precisely, the half plane p = with H > is a repeller, whereas 
when H > it is an attractor. This means that, at late time, all the solutions go towards this half plane. Moreover, in 
the contracting phase there is a critical point ( — y/ — l/f3, 0, 0) (the contracting de Sitter universe) which restricted to the 
plane p — is a repeller. This means that only a solution tends asymptotically towards this point, while all the other 
escape towards infinity in finite time (this was proven in [ 1 1). The same does not happen for a phantom field where the 
manifold p = is not invariant, and eventually, the system crosses this manifold, i.e., it can have negative energy density, 
and then it cannot leave the decreasing phase, becoming singular at finite time, as we have shown, again numerically and 
analytically. 



9 



Acknowledgments. This investigation has been supported in part by MICINN (Spain), projects MTM201 1-27739-C04- 
01, MTM2009-14163-C02-02, and FIS2010-15640, by the CPAN Consolider Ingenio Project, and by AGAUR (Gen- 
eralitat de Catalunya), contracts 2009SGR 345, 994 and 1284. EE was also supported by MICINN (Spain), contract 
PR201 1-0128, and his research was partly carried out while on leave at the Department of Physics and Astronomy, Dart- 
mouth College, 6127 Wilder Laboratory, Hanover, NH 03755, USA. 



[1] J. Haro, J. Amoros and E. Elizalde, Phys. Rev. D83, 123528 (2011). 

[2] S. Perlmutter et al., Astrophys. J. 517, 565 (1999). 

[3] A.G. Riess et al., Astron. J. 116, 1009 (1999). 

[4] E. Komatsu et al., Astrophys. J. Suppl. Ser. 192, 18 (2011). 

[5] S.K. Srivastava, Gen. Relativ. Gravit. 39, 241 (2007). 

[6] S. Nojiri, S. Odintsov and S. Tsujikawa, Phys. Rev. D71, 063005 (2005). 

[7] V. Mukhanov, Physical fundations of cosmology (Cambridge University Press, 2005). 

[8] M.V. Fischetti, J.B. Hartle and B.L. Hu, Phys. Rev. D20, 1757 (1979). 

[9] P.C.W. Davies, Phys. Lett. B68, 402 (1977). 
[10] S. Nojiri and S. Odintsov, Phys. Rev. D70, 103522 (2004). 
[11] E. Elizalde, S. Nojiri and S. Odintsov, Phys. Rev. D70, 043539 (2004). 
[12] T. Azuma and S. Wada, Prog. Theor. Phys. 75, 845 (1986). 
[13] S. Wada, Phys. Rev. D31, 2470 (1985). 

[14] A. Ashtekar and P. Singh, Class. Quantum Grav. 29, 213001 (201 1). 

[15] P. Singh, Class. Quantum Grav. 26, 125005 (2009). 

[16] P. Singh, J. Phys. Conf. Ser. 140, 012005 (2009). 

[17] T. Naskar and J. Ward, Phys. Rev. D76, 063514 (2007). 

[18] D. Smart and B. Gumjudpai, Phys. Rev. D76, 043514 (2007). 

[19] G. Cognola, E. Elizalde, S. Nojiri, S. Odintsov and S. Zerbini, Phys. Rev. D73, 084007 (2006). 

[20] J. de Haro, Future singularity avoidance in phantom dark energy models (to appear in JCAP), gr-qc: 1204.5604 (2012). 



